Impurity spin relaxation in S = 1/2 XX chains 
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Dynamic autocorrelations {S°'{t)S") {a = x,z) of an isolated impurity spin in a 5* = 1/2 XX 
chain are calculated. The impurity spin, defined by a local change in the nearest-neighbor coupling, 
is either in the bulk or at the boundary of the open-ended chain. The exact numerical calculation 
of the correlations employs the Jordan-Wigner mapping from spin operators to Fermi operators; 
effects of finite system size can be eliminated. Two distinct temperature regimes are observed in 
the long-time asymptotic behavior of the bulk correlations. At T = only power laws are present. 
At high T the x correlation decays exponentially (except at short times) while the z correlation still 
shows an asymptotic power law (different from the one at T = 0) after an intermediate exponential 
phase. The boundary impurity correlations ultimately follow power laws at all T, with intermediate 
exponential phases at high T. The power laws for the z correlation and the boundary correlations 
can be derived from the impurity-induced changes in the properties of the Jordan-Wigner fermion 
states. 
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I. INTRODUCTION 



The S =1/2 XX chainl 



Af-l 



(1.1) 



is one of the simplest quantum many-body systems con- 
ceivable, as many of its properties can be derived from 
those of noninteracting lattice fermions. Its equilibrium 
spin pair correlation functions 



Tre^^^ 



(1.2) 



have been the objects of intense research efforts over an 
extended periodcl Ell. Only a few explicit analytic results 
are available, but several existing asymptotic results for 
large distances |« — j| or long times t have been corrobo- 
rated by numerical calculations. 

For XX chains with homogeneous nearest-neighbor 
coupling (Ji = J) only three different types of asymp- 
totic long-time behavior have been observed to date: 
Gaussian, exponential, and power-law (often with su- 
perimposed oscillations). It is interesting to speculate 
whether non-uniform or random couplings might induce 
additional types of asymptotic behavior. 

In the present paper we study the changes in autocor- 
relation functions {S°'{t)S°') induced by a single impurity 
spin in an otherwise homogeneous chain. The impurity 
spin is located either at the boundary of the system, 



Ji = J', J, = J = 1 for i > 2 



(1.3) 



or in the bulk, 



J, 



N/2-1 



J, 



N/2 



J', = J 1 for all other i. (1.4) 



Equilibrium and non-equilibrium dynamics of, the 
boundary impurity were studied early on by TjonB. In 



the weak-coupling limit (J' — > 0), Tjon obtained ex- 
ponential behavior of the impurity spin autocorrelation 
functions. Our results (see Sec. IV) show that for finite 
impurity coupling J', exponential behavior occurs only in 
an intermediate time regime, whereas the ultimate long- 
time behavior is a power law. 

Besides their obvious relevance in low-dimensional 
magnetism, impurities in spin-1/2 chains are also of in- 
terest in quantum dynamics, where two- level systems 
coupled to "baths" serve as rriadals for quantum sys- 
tems iiLjdissipative environmentsc3~cil. The most popular 
modelcj in this field is the spin-boson model, consisting 
of a single spin 1/2 coupled to a (quasi-) continuum of 
noninteracting oscillators with a given spectral density. 
In a recent studynJ the oscillator bath was replaced with 
a bath of noninteracting spins 1/2. The changes in dy- 
namic behavior which were observed as a result of this 
replacement suggest further exploration of different kinds 
of baths. The system studied here can be considered a 
two-level system (the impurity spin) coupled to a bath of 
interacting two-level systems (the remainder of the XX 
chain). An interesting feature of this system is the fact 
that while the z component of the total spin is conserved, 
the X component is not. Thus differences are to be ex- 
pected between the relaxation of the x and z components 
of the impurity spin. 

The plan of the paper is as follows: In Sec. || we dis- 
cuss the method used to calculate the dynamic correla- 
tion functions numerically. In Sec. Ill we present results 



for spin autocorrelation functions of a bulk impurity spin 
(and also of its neighbors) for both zero and finite T. In 
Sec. we discuss boundary impurity autocorrelation 
functions. Sec. |^ summarizes our findings. 
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II. METHOD 

The open-ended iV-site spin-1/2 XX chain described 
by the Hamiltonian (^]^) can be mapped to a Hamilto- 
nian of noninteracting fermions, 



N-l 



(2.1) 



Ek (2.4) with J'-dependent k values. For J' > Jc a, pair 
of exponentiaUy locahzed impurity states with energies 
±eoi ko| > Ij emerge from the continuum. The critical 
coupling strength is Jc = 1 for the bulk impurity and 
Jc — \/2 for the boundary impurity. Below, we shall oc- 
casionally refer to properties of the analytic solution in 
order to explain the long-time asymptotic behavior ob- 
served in the numerical results. 



i=l 



by means of the Jordan- Wigner transformationflfl be- 
tween spin and fermion operators: 



1 



(2.2) 



III. BULK IMPURITY 



A. T = 



S+ = (-1)E;:.4=.4 = _ 2clck)cl (2.3) 



fe=l 



In the homogeneous case Ji = J, the one-particle energy 
eigenvalues are 



Ek — — J COS k,k = 



N + 1 



(2.4) 



and the eigenvectors are sinusoidal functions of the site 
index i. For general Ji neither eigenvalues nor eigenvec- 
tors are available analytically, however, both are easily 
obtained from the solution of a tridiagonaLeigenvalue 
problem with standard numerical proceduresEEI. 

The spin correlation functions (1.2) are mapped to 
fermion correlation functions, with crucial differences be- 
tween the cases a = z and a — x. {Sf{t)Sj) maps to a 
density-density correlation function involving four Fermi 
operators. Due to the string of signs in (2.3), however, 
{Sf{t)Sj) maps to a many-particle correlation function 
involving 2{i + j — 1) Fermi operators. Wick's theorem 
can be applied to expand {Sf{t)Sj) in products of el- 
ementary fermion expectation values. That esBansion 
can be most compactly expressed as a Pfaffiancll whose 
elements are sums involving the one-particle eigenvalues 
and eigenvectors. 

In order to obtain results valid in the thermodynamic 
limit N oo, finite-size effects must bejidentified and 
eliminated. As finite-size effects are knownEiJ to be caused 
by reflections of propagating excitations from the bound- 
aries of the system, the maximum fermion group velocity 
(see 2.4) can be used to estimate the time range over 
which a given spin correlation function (1.2) can be ex- 
pected to be free of finite-size effects. That estimate 
can then be verified by explicit numerical calculation of 
equivalent correlation functions for system sizes iVo and, 
say, 2No. 

The one- fermion eigenvalue-, problem for the single- 
impurity chain ( |1.3D or (1.4)c3 may be solved analyti- 
cally. The nature of the solution depends on the value 
of J'. For J' below a critical value Jc all states are ex- 
tended and the continuous energy spectrum is given by 
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FIG. 1. Bulk impurity spin x autocorrelation function 
l(*S'^/2(*)5'^/2>l at T = 0. Main plot: N = 512, impurity 
coupling J' < 1. The heavy solid lines (top to bottom) cor- 
respond to J' = 0.1, 0.15, 0.2, 0.4, O^and 0.8. The dashed 
straight line is the t~^^^ power law (3.1). The thin solid lines 
are A'^ — 256 data demonstrating the influence of system size. 
Inset: A'^ = 128, J' > 1. Shown are data for J' = 1 (heavy 
solid line), J' — 1.2 (Dashed curve), and J' — 1.4 (thin solid 
I me) . The dashed straight line is again the t~^^^ power law 

The long-time asymptotic behavior of the T ~ bulk 
impurity spin x autocorrelation function is difficult to 
obtain due to a combination of two reasons. Firstly this 
correlation is the computationally most demanding one, 
as large PfafRans have to be evaluated. Secondly it is also 
the correlation function displaying finite size effects at the 
earliest times. This may be related tp . its particularly 
slow long-time asymptotic decay lawOiiJIlj 



(^f W^.^n) - in' Jh^)-"^ for T = (3.1) 



2 



in the homogeneous case Ji = J. It should be noted 



that the right-hand side of ( 3J ) is the leading term of an 
asymptotic expansion; its character changes from purely 
real (for J'^t^ < -n?) to complex (for J'^t^ > n^). More 
explicit forms are eq. (1-23) in ref. ^ and eqs. (59,61) 
in ref. |l^. 

We have calculated {S%/2(^)^n/2) impurity cou- 
pling constants 0.1 < J' < 4. In all cases the asymptotic 
decay of th e co rrelation function was consistent with the 
law With growing J' \{S-^/^{t)S%^2)\ de- 

velops oscillations of rather well-defined frequency and 
growing amplitude, as shown in the inset of Fig. ^ The 
frequency of the oscillations for J' > 1 is proportional to 
the energy 



£0 



J'- 



12 



1 



(J' > 1) 



(3.2) 



of the localized impurity state. On the whole, the long- 
time asymptotic behavior of the impurity spin x auto- 
correlation at T = is not fundamentally changed by 
varying the value of J'. 




FIG. 2. Main plot:BuIk impurity spin z autocorrelation 
function K5'^/2(t)S'^/2)l {N = 1024) at T = 0, for impu- 
rity couplings J' — 1,0.99,0.8 (solid lines, top to bottom); 
J' = 0.6, 0.4, 0.2 (long-dashed lines, bottom to top). The thin 
straight lines show the power laws and . Inset: same 
as main plot, for J' = 1, 1.2, 1.4, 2, 4 (bottom to top). 

The impurity spin z autocorrelation, in contrast, 
changes significantly when J' is varied, as shown in Fig. 
^. For small J' \{S^^2i^)^N/2)\ displays a monotonic de- 
cay. Roughly at J' = 0.5 oscillations (of well-defined 
and J'-independent frequency) begin to develop. For all 
J' < 1 the correlation function follows an asymptotic 



la-aj—IEhe J' = 1 correlation shows the law well 
knowrJj'ErEj for the homogeneous case. As soon as J' 
is further increased, the behavior changes again and the 
absolute value of the correlation function tends to a con- 
stant nonzero value for large times. 

The changes in the asymptotics of |(<S'^/2(^)'^Af/2)l '^^^^ 
be understood from the analytic solution mentioned in 
Sec. ||. The Jordan- Wi gner transformation yields, after 
a few simple steps, 

{S!{t)S!)=(^m:^)\'e^'^'f{e,)^ (3.3) 



Here |j/) is a one-fermion eigenstate of H ( PT| ) with en- 
ergy El, and f{x) = (exp(/3a;) -1-1)"^ is the Fermi func- 
tion. For J' > Jc — 1 the presence of a localized im- 
purity state with large and with outside the 
continuum yields a harmonically oscillating non-decaying 
contribution to (S'f (t)S'f ). Similar contributions are con- 
tained in every element of the Pfafhan for {Sf{t)Sf), but 
not in that correlation itself (see Fig. |^, inset). The 
reason probably is a cancellation of terms due to the 
multiplications and additions inherent in the definition 
of the Pfafhan. T he t ime scale introduced by the dis- 
crete energy value (|3.2|j is reflected in the oscillations of 
K'^Af/2(^)'^Af/2)l (Fig. y inset). Similar behavior is found 
for T > and will be discussed in the next subsection. 

For J' < I the time-dependent term in ( 0.3| ) is propor- 
tional to 



de(l-e^)-^e-7(e)l(*|e)l 



(3.4) 



where {i\e) corresponds to (i|i^) in (3.3) and the inverse 
square root factor is the one-particle density of states 
of the dispersion (^^) (which still describes the energy 
eigenvalues, only with slightly displaced k values for J' ^ 
1). For J' — 1, |(i|e)P docs not depend on e, the inverse 
square root singularities at the band edges e ^ ±1 lead 
to a t~^/^ asymptotic behavior of the integral, and to a 
t-i behavior of {Sf(t)Sf). For J' < 1 the amplitude of 
the one-particle eigenstate with energy e at the impurity 
site is 



J'- 



1- J'- 



J'2 



(3.5) 



(apart from weakly e-dependent normalizatio n fa ctors). 
This changes the band-edge singularities in (^^) from 
(1 — e^)^^/^ to (1 — e^)^/^, so that the integral contains a 
t^'^/^ term. Consequently, {Sf{t)Sl) contains &t~'^ term 
which dominates for T > (see next subsection). At 
T = 0, however, the leading term is proportional to t ~'^ 
due to the discontinuity of the Fermi function in (3.4). 
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We have also studied the x and z autocorrelations of 
nearest and next-nearest neighbors of the impurity spin 
i = N/2. For weak impurity coupling {J' < 0.2) the x 
correlation functions of spins i — N /2 + 1 and i — N/2 + 2 
show weak oscillations superimposed on a t^^^^ decay 
masked by strong finite-size effects. The z correlations 
(for J' < 1) show stronger oscillations. Their decay looks 
roughly like a power law with an exponent between -2 and 
-3. 



first and later crosses over to the law derived above, 
with superimposed oscillations which are absent in the 
low-temperature regime. The appearance of oscillations 
(if only of small amplitude) in (S'f (i)S'f) at hifih T is 
reminiscent of the phenomena recently reportcdcS for a 
two-level system coupled to a spin bath. In that sys- 
tem, a persistence of oscillations up to infinite T could 
be observed. 



B. T > 




200 



FIG. 3. Crossover between the regimes of high tempera- 
ture (lower set of curves) and low temperature (upper set of 
curves). Main plot: Impurity spin autocorrelation function 
\{SN/2i't)SN/2)\ (N = 256) with impurity coupling J' = 0.3 
for temperatures T = 10^ 10*, 10^ 10°, 10~\ 10"^ 10"^ 
10"", 10"^ and 10"^ (bottom to top). The curves for 
T = 10"^ and 10"^ are long-dashed, all others are solid. Note 
that several curves coincide on the scale of the figure. Inset: 
same as main plot, for |(<S'^/2(^)'S'Ar/2)l i-^ ~ 128); time range 
is < t < 100. 

Whereas at T ^0 only power-law decay is observed, 
exponential decaja becomes possible at finite T. Fig. ^ 
shows X and z impurity spin autocorrelations at J' — 0.3 
for several decades in T. There are two well defined tem- 
perature regimes with a crossover between them. Within 
each regime the correlation functions do not change qual- 
itatively: note that several of the curves in Fig. ^ coin- 
cide. The X autocorrelation in the high-T regime shows 
exponential decay which persists over the entire time 
range during which the the results are free of finite-size 
effects. The z autocorrelation decays exponentially at 
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FIG. 4. Exponential decay of the bulk impurity spin x 
autocorrelation function at high temperature. Main plot: 
|(5'^/2(i)S'^/2)l = 128) at T = 10^ for impurity couplings 
J' = 0.1,0.15,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9, and 1 (top to 
bottom). Inset: Exponential decay rates determined from 
the data shown in the main plot, and from analogous data 
for J' > 1 and for T = 1. The circles are T = 10^ data, the 
crosses are T = 1 data. Note the pronounced differences near 
J' — 1, where the infinite-T decay rate diverges. 

Fig. I shows |(5'^/2(i)5'^/2)l at T = 10^ for impurity 
coupling 0.1 < J' < 1 in a TV = 128 chain. As autocor- 
relations at T = oo are real even power series in t, all 
curves start with zero slope at t = 0, but then (with the 
exception of J' = 1) bend over to a nearly perfect ex- 
ponential decay. The inset shows the decay rate of that 
exponential decay as fitted to the data in the main plot. 
Also shown is the decay rate determined from data for 
T = 1, and for J' > 1. Differences between T — 1 and 
T — oo are to be expected and are indeed visible in the 
behavior of the decay rate close to J' = 1: For finite T 
the X autocorrolation function of a homogeneous chain 
Ji = 1 is knownc^Eil to decay exponentially with a finite 
T-dependerit|-decay rate, whereas for infinite T the decay 
is GaussiarJacl^tSl. The Gaussian decay corresponds to a 
divergence of the exponential decay rate which is obvious 
in the inset of Fig. |4[ 
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For J' > 1 \ {S^i2{t)S^i^ \ is no longer (almost) purely 
exponential, but develops considerable oscillations. The 
exponential decay rate grows with T, but decreases as J' 
grows, as shown in the inset of Fig. |^. As in the T = 
case the frequency of the oscillations is proportional to 

In order to obtain a quantitative measure of the preci- 
sion to which the decay of \{S'^ i^{t)S'^ at T = 10^ 
follows an exponential, we fitted an exponential law 
aexp(— &/:) to the numerical data for < i < 100 and 
calculated the quantity 



pit) 



(3.6) 



a cxp{—bt) 

which should equal unity for a purely exponential decay. 
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FIG. 5. Precision of the exponential decay of the 
high-temperature bulk impurity spin cor relation function. 
Shown is the precision function p{t) ( |3.6| ) for the data for 
J' < 0.4 from the main plot of Fig. |. 

For J' < 0.4, p{t) is shown in Fig. ||. Note that the 
scale of the figure extends only to a maximum deviation 
of 1 percent from purely exponential decay. The general 
slope in the data is a natural consequence of the intrin- 
sically non-exponential behavior of the correlation func- 
tion for small t which was mentioned above. The data 
for J' < 0.2 follow the purely exponential fit to a preci- 
sion of better than 2 parts in thousand for t > 10. This, 
rules out the stretched-exponential behavior reportedE2l 
for J' < 0.2 at T = oo in an approximate study based 
on extrapolation of truncated continued-fraction expan- 
sions. Fig. H also reveals the presence of tiny oscillations 
which are invisible on the scale of Fig. ^. The frequency 
of these oscillations is independent of J', in contrast to 



the stronger oscillations for J' > 1 already mentioned 
above. 
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FIG. 6. Initially exponential and asymptotically power-law 
decay of the bulk impurity z spin autocorrelation function at 
high temperature. Main plot: |(S'^/2(i)S^/2)l (N = 256) 
at T = lO'^ for impurity couplings J' = 0.15,0.2, and 0.3 
(top to bottom). Note how the exponential regime gets 
shorter as J' grows. The features in the lower right corner 
are finite-size effects. Inset: Maxima of |(>S'^/2(i)'S'^/2)l for 
J' = 1, 0.9, 0.8, 0.7, and 0.6 (top to bottom) and T > 20. The 
dashed straight lines represent the and power laws. 

The main differences between x and z spin autocor- 
relation functions were already shown in Fig. ^ Some 
more detail on the behavior of |(<S'^/2(0'5'Ar/2) I P''^" 
sented in Fig. ^. The main plot shows the crossover 
between the exponential and (with superimposed os- 
cillations) regimes for three small J' values. With grow- 
ing J' the exponential decay rate grows and the expo- 
nential regime shortens in such a way that in the expo- 
nential regime |(S'^/2(^)'^Ar/2)l ^ decreasing function of 
J' whereas in the power-law regime it is an increasing 
function of J'. We have deliberately chosen a very long 
time window in order to illustrate how finite-size effects 
manifest themselves (for t > 230). In the inset of Fig. ^ 
we demonstrate the asymptotic power-law behavior for 
larger values of J'. The dots represent the maxima of 
\{S^N/2it)SN/2)\^ which follow the t'^ (for J' < 1) and 
(for J' = 1) laws already discussed. For J' > 1 the 
behavior changes to a "constant with oscillations" type 
of asymptotics, similar to the T = situation shown in 
the inset of Fig. |[ For T > 0, however, the amplitude 
of the oscillations is considerably larger than for T = 0. 

The high-temperature spin autocorrelations of the 
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nearest and next-nearest neighbors to the impurity (for 
J' < 1) do not show any particularly surprising features. 
The z correlations do not show exponential decay in the 
beginning. Instead they oscillate and the maxima of the 
oscillations display the familiar (for J' < 1) and t^^ 
(for J' = 1) laws. The x autocorrelations interpolate 
smoothly between two known limiting cases. At J' = 1 
the X autocorrelation of the spin i = N/2 + 1 of course 
is a Gaussian as that of any other bulk spin. At J' = 0, 
however, i = N /2-\-l is the first spin in a semi- infinite ho- 
rrpgeneous chain, whose T — oo autocorrelation function 
ista a combination of Bessel functions with an asymp- 
totic t~^/^ decay. Upon reducing J' from 1 to 0, the 
development of the characteristic Bessel function oscilla- 
tions (with zeros hardly depending on J') can be nicely 
observed. Similarly the time range during which the cor- 
relation function follows the expected t'^^"^ decay grows 
as J' diminishes. The x autocorrelation of ? = N/2 + 2 
behaves quite similarly, but only a small number of os- 
cillations is visible (in a linear plot) because of the fast 
^^-9/2 -J asymptotic decay of the knowro J' — Bessel 
function expression. 



IV. BOUNDARY IMPURITY 



The boundary impurity is defined by (1.3). Similarly 
to the case of the bulk impurity discussed in the previous 
section, the boundary correlation functions show a low- 
temperature regime and a high-temperature regime, and 
we restrict our discussion to the values T — and T = 1 
which represent these two regimes. 

It suffices to discuss the impurity spin x autocorrela- 
tion function {Sf{t)Sf), because 



2(5f(05r) = 5]i(ii 



(4.1) 



that is, the square root of the time-dependent part of 
{Sf{t)Sf) (see O). The presence of an isolated impu- 
rity state for J' > = \/2 (compare Sec. |l|) should 
be visible in the dynamic correlation functions. In fact, 
the asymptotic behavior of {Sf{t)Sf) displays a crossover 
similar to the one shown in the inset of Fig. 13: for all T 
the long-time behavior is a power law for J < \f2 and 
a constant for J' > v^, with additional oscillations in 
both regimes. 
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FIG. 7. Boundary spin x autocorrelation function 
K^i (t)S'i )|. Main plot; High-temperature regime, T = 1, 
iV = 512, and impurity coupling J' = 1, 0.8, 0.6, 0.4, 0.2, 
and 0.1 (left to right, according to the point where the curves 
enter the figure at the upper boundary). The dashed line 
shows the t''^^"^ law. Inset: Low-temperature regime, T = 0, 
iV = 1024, and J' = 1, 0.7, 0.5, and 0.1 (bottom to top). The 
dashed line shows the law. 

Fig. 1^ shows |(S'f (i)S'f )| for impurity coupling J' < 1. 
For r = 1 we observe an initially exponential decay fol- 
lowed by a power law. The exponential decay rate grows 
with J' . The duration of the initial exponential phase 
decreases with growing J' in such a way that in the sub- 
sequent power-law regime the correlation function is an 
increasing function of J'. This behavior is (not unex- 
pectedly, compare (4.1)) similar to that of the bulk im- 
purity z autocorrelation discussed in Sec. |IIIB| , however, 
the asymptotic power law is a different-one, namely the 
t~'^/^ law known for the boundary spinll3 of the homoge- 
neous semi-infimte chain at infinite T, and for a range of 
boundary spinsEll of the same system at finite T . _ 

At this point, the early analytical study by Tjontl 
should be mentioned. In the limit of sufficiently small J' 
Tjon found an asymptotically exponential decay ^ e^*/'^ 
for (S'f (t)S'f), with decay rate = ^J'^. Indeed, our 
numerical data show that during the exponential regime 
mentioned above, the decay rate is quite precisely equal 
to I J'^ for J' < 0.2, and a bit larger for larger J'. How- 
ever, we also observe numerically (and explain analyti- 
cally, see below) a crossover from exponential to power- 
law behavior. The duration of the exponential regime 
grows (~ J'^^) as J' becomes weak and thus our numer- 
ical results for arbitrary J' connect smoothly to Tjon's 
analytical result restricted to the weak-coupling limit. 

The inset of Fig. @ shows \{Sf{t)Sf)\ for T = and 
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J' < 1 . After a very slow initial decay (slower for smaller 
J') the curves eventually all bend over to show a, 
asymptotic decay developing some oscillations as J' ap- 
proaches unitji—The decay at T = is again a well- 
known featurdl2l of the homogeneous semi-infinite chain. 

The asymptotic power laws may be understood from 
the properties of the analytic solution of the boundary 
impurity (one-par ticle) problem, along the lines of the 



seen in further studies. 



discussion in Sec. |IIIA| . According to (yj), \{Sf{t)Sf)\ 
for J' <V2 is given by an integral analogous to the one in 
( |3.4[ ). The analytic solution shows that the wave function 
factor in the integral, (Ije) ~ sinfc, where e = — cosfc. 
This leads to a band-edge singularity ~ (1 — e)^/^ in 
the integrand, and a t~^^'^ asymptotic behavior of the 
integral. At T = the contribution is dominated by 
the t^^ contribution from the discontinuity of the Fermi 
function. 





low T 


high T 


X 


^-1/2 


t 

e -r 


Z,J' <1 










(after initial e~~) 


Z,J' >1 







TABLE I. Asymptotic decay of the bulk x and z autocor- 
relation functions. Additional oscillations of varying strength 
are present in all cases. 



V. SUMMARY AND CONCLUSIONS 



The dynamic spin correlation functions associated with 
isolated impurities in a 5 = 1/2 XX chain show a rich be- 
havior depending on the temperature, the impurity cou- 
pling strength J', the spin component (a = x, z) under 
consideration, and on the position of the impurity spin 
in the chain. Regimes of low and high T, with qualita- 
tively different behavior of the correlations may be distin- 
guished. We have summarized the asymptotic behavior 
of the bulk x and z spin autocorrelations in these two 
T regimes in Table |. The basic features of the bound- 
ary correlations may also be obtained from Table | by 
observing that (i) the z autocorrelation does not change 
fundamentally between bulk and boundary and (ii) at 
the boundary the x autocorrelation behaves as the square 
root of the z autocorrelation. 

As mentioned in the introduction, the present model 
may be interpreted as a single two-level system in con- 
tact with a large bath. The influence of the nature of the 
bath degrees of freedom on the typie^of decay of the two- 
level system was addressed recentlycj for a spin bath, con- 
structed in a way to resemble closely the standards har- 
monic oscillator bath. The following results were found. 
At T = the spin bath leads to damped oscillations in 
the two-level system, as does the oscillator bath. How- 
ever, at high r, the oscillations vanish for the oscillator 
bath but persist for the spin bath. 

Without going into any detailed comparison between 
the spin bath employed in Ref. ^ and the XX chain 
studied here, we would like to point out the existence of 
similar oscillation phenomena at high T in the present 
system. Fig. || (main plot) shows how the bulk impu- 
rity spin z autocorrelation function develops oscillations 
as temperature increases. Whether these oscillations can 
be unambiguously assigned to either the impurity or the 
bath, and what happens for systems interpolating be- 
tween the present one and that of Ref. Eq, remains to be 
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